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SUMMARY 


Several approximate procedures for calculating the bending -moment 
response of flexible airplanes to continuous isotropic turbulence are 
presented and evaluated. The modal methods (the mode -displacement and 
force -summation methods) and a matrix method (segmented-wing method) are 
considered. These approximate procedures are applied to a simplified 
airplane for which an exact solution to the equation of motion can be 
obtained. The simplified airplane consists of a uniform beam with a 
concentrated fuselage mass at the center. Airplane motions are limited 
to vertical rigid-body translation and symmetrical wing bending deflec- 
tions. Output power spectra of wing bending moments based on the exact 
transfer -function solutions are used as a basis for the evaluation of 
the approx ima te methods. It is shown that the force -summation and the 
matrix methods give satisfactory accuracy and that the mode -displacement 
method gives unsatisfactory accuracy. 


INTRODUCTION 


Current analytic methods of studying airplane responses to con- 
tinuous atmospheric turbulence make use of the techniques from random- 
process theory. These techniques provide a compact description of the 
statistical properties of the turbulent velocities in a form suitable 
for response calculations. The significant characteristics of turbu- 
lence velocities for this purpose are contained in power spectra and 
appropriate cross-spectra. These spectra are then used along with the 
appropriate airplane frequency-response or transfer functions for deter- 
mination of the airplane response. 

The early applications of this technique described the turbulence 
as a one -dimensional random process varying along the flight path of the 
airplane but uniform along the wing span. (See, for example, ref. 1.) 

In recent studies the variation of turbulence across the wing span 
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(two-dimensional or isotropic turbulence) has been shown to give rise to 
significant effects for large flexible airplanes. (See, for example, 
ref. 2.) 

In general, the solution of the dynam:.c equations of motion for the 
transfer functions requires the use of approximate procedures. It is 
the purpose of the present report to investigate the accuracy of several 
of these approximate procedures used for calculating the bending -moment 
transfer functions when applied to the case of a flexible airplane sub- 
jected to two-dimensional turbulence. The approximate methods under 
consideration are the mode -displacement method, the force -summation 
method, and a matrix method ( segmented -w in*; method - also referred to as 
lumped-parameter method) . An account of these approximate procedures 
can be found in reference 3, which also cites other references. 

The approximate procedures under investigation are applied to a 
simplified airplane for which an exact solution to the equation of motion 
can be obtained. These exact-solution results serve as a basis for the 
evaluation of the approximate methods. For purposes of comparison, 
results are also calculated for the f lexib Le airplane traversing one- 
dimensional turbulence and a rigid airplane traversing one- and two- 
dimensional turbulence. 

The simplified airplane consists of a uniform beam with a concen- 
trated fuselage mass at the center. Airplane motions are limited to 
vertical rigid-body translation and symmetrical wing bending deflections 
only. Aerodynamic forces are based on strip theory and include the 
effects of unsteady lift associated with tie gusts and the resulting 
wing motions. 

A secondary purpose of this report is to indicate some effects of 
spanwise variations of turbulence in addition to those given in refer- 
ence 4. This purpose is accomplished by comparing mean-square bending- 
moment responses of various airplane configurations to one- and two- 
dimensional turbulence. The airplane configurations are defined by the 
ratio of fuselage mass to total mass. 


SYMBOLS 


b wing span, ft 

C(k) generalized Theodorsen function 

wing lift-curve slope per radian 
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c local wing chord, ft 

c average chord, ft 

c 7 section lift-curve slope per radian 

El bending stiffness of wing, lb/sq ft 

F(k), G(k) real and imaginary parts of generalized Theodorsen function 
C(k), respective ly 

H(k) transfer function of bending moment for one-dimensional 

turbulence 


H(k, q) , H( k, q+9) influence transfer functions of bending moment for 

two-dimensional turbulence 


Hq transfer function of bending moment resulting from vertical 

wing motion in rigid-body mode 

transfer function of bending moment resulting from vertical 
wing motion in fundamental wing bending mode 

H 2 transfer function of bending moment resulting from external 

gust loading 

K 0 ( ), K l ( ) Bessel functions of second kind, order 0 and 1 

CJJC 


reduced frequency, 


2U 


reduced frequency associated with fundamental wing bending 

(i)|C 

mode , 

’ 2U 


k a = 


(c/ 2 ) 2 EI 


(b/2) 5 U 2 (m w /2) 


L 

lc 


scale of turbulence, ft 
gust load per unit span 


inertia load per unit span 
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“m 


B 


M 
M B,s 

m(y) 

m p 

my 


aerodynamic load per unit span due to wing vertical motion 
bending moment, ft -lb 

bending moment due to sharp-edge gust, — Ct qSy CD — , ft-lb 

2 a F n 

wing mass distribution, slugs/ft 
fuselage mass, slugs 

wing mass, slugs 


m^(y), m^.(y) total airplane mass distributions, slugs/ft 


m t 

q 

R 

Re( ) 

S 

s 

t 

U 

w 

y 

y 

^cp 


- - y °P 

Y °P b/2 

yi 


“f _ H 

mw 1 - M. 


total airplane mass, slugs 
dynamic pressure, lb/sq ft 
ratio of fuselage mass to wing mass, 

designates real part of ( ) 
wing area, cb, sq ft 

dimensionless scale of chordwise turbulence. — 

c/2 

time, sec 

forward velocity, fps 

vertical component of gust velocity, fps 
spanwise displacement, ft 
dimensionless spanwise displacement, 

b/1 77 


lateral ordinate of center of pressure, - y^ + — yp j, ft 


1 77 2 


spanwise location of bending -moment station, ft 
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yi 


z 


z 


r(e) 

7 n( ) 
8 ( ) 
’ll 


9 


K 


n 

6 0 < ) 
Sl< > 

P 

a 2 

®w( k ) 


dimensionless spanwise location of bending -moment station, 

yi 


b/2 


vertical displacement, ft 

dimensionless vertical displacement, — ^ — 

c/2 

gamma function (see eq. (13)) 
space function used in equation (9) 


Dirac delta function 
spanwise coordinates, ft 

dimensionless spanwise coordinate, — — 

b/2 


variable of integration, - T), ft 

8mt 


airplane mass parameter, 


C lu pSc 


wing mass parameter, 


fuselage mass ratio, — 

m t 


C^pSc 

mg 


= k( 1 - ji) 


rigid-body mode 
fundamental wing bending mode 


air density, slugs/cu ft 

mean-square deviation of bending moment, lb -ft 
output power spectrum of bending moment 

power spectrum of vertical gust velocity for one -dimensional 
turbulence 
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$ w (k,e), $ w (k, t^-tA power spectra of vertical gust velocity for tvo- 
' ' dimensional turbulence 


(k) 

w mn 


Subscripts: 


Notation: 


' dimensional turbulence 
associated gust spectra (see cq. (12)) 

unsteady lift function for gust penetration (Sears function) 
frequency, radians/sec 

frequency associated with fundamental wing bending mode, 
radians/ sec 

unit step function 


flexible airplane 

imaginary component of complex variable 

integers, 0, 1, 2, . . . 

real component of complex variable 

Rayleigh 

rigid airplane 


* complex conjugate 

| | absolute value 

[ 1 square matrix 

[* J diagonal matrix 

[ J row matrix 

{ } column matrix 

A dot over a symbol indicates a derivative with respect to time, 
and a prime indicates a derivative with respect to y. 
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BASIC INPUT-OUTPUT RELATIONS 


Random-process techniques as applied to the gust-load problem are 
discussed in many references. (For example, see refs. 1 to 4.) From 
reference 4, the output spectra of bending moments are defined as follows 

For one -dimensional turbulence, 

*MB(k) = | H(k)| 2 1> w (k) (1) 


where 


output power spectrum of bending moment 


H(k) transfer function for bending -moment response of airplane to 

Continuous sinusoidal gusts of unit amplitude and wave 
length jtc/k 

<D (k) input power spectrum of vertical gust velocity 


For two-dimensional turbulence, 


r b/2 f b/2 

*MB(k) = / / H(k,Tj)H*(k,ii 1 )« w (k,Ti 1 -Ti)dTi din (2) 

J_b/2^-b/2 


where 


H(k, rj) influence transfer function which represents the complex 

amplitude of bending -moment response to continuous sinusoi- 
dal gust of unit amplitude and wave length nc/k impinging 
at station rj (* indicates complex conjugate) 




input power spectrum of vertical gust velocity 


It has been shown in reference 4 that equation ( 2 ) may be simplified 
somewhat by a change of variable; namely, = 0 + t]. Thus, performing 

the integrations first over tj and then over 9, as shown in reference 4 
and restricting the system to the symmetrical response only yields 



(3) 
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where Hp( ) and Hj( ) are, respectively, the real and imaginary 
components of the symmetrical transfer functions. 

The output spectra, defined by equations (l) and ( 3 ), can in turn 
be used to calculate many other statistical quantities (see ref, 5); 
however, the only other statistical quantity to be calculated herein is 
the mean-square deviation of the bending moment. From reference 4, the 
mean-square response is defined as 

o 2 = f %B( k ) dk: (*0 

J 0 


TRANSFER FUNCTIONS 


From simple beam theory, the basic differential equation of bending 
motion is 


a 2 

dy 2 




U + ^ir. + 


(5) 


where l ^ is the inertia load per unit span and l g and Z m are, 

respectively, the aerodynamic loads per unit span due to the disturbing 
gust and those due to wing vertical motions. 

Solutions of equation (5) to obtain airplane motions and deflections 
and, subsequently, the determination of the bending-moment transfer func- 
tions are accomplished by an exact method, two modal methods (mode dis- 
placement and force summation), and a matrix method (segmented-wing 
method) . 


All of the solutions are made for the simplified airplane, consisting 
of a uniform beam with a concentrated fuselage mass at the center, shown 
in figure 1. It is assumed that the airplane is restricted to vertical 
motions only and, in addition, that the pitching motions be neglected. 


Exact Method 

Equation ( 5 ) is solved by means of the Laplace transform method to 
obtain z M and, thence, the bending-moment influence transfer function 
from EIz" . A brief derivation of the exact analysis is given in 
appendix A for one- and two-dimensional turbulence. 
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Modal Methods 

The fundamental assumption in all modal methods is that the displace 
ment can be approximated by the superposition of a finite number of nat- 
ural or assumed mode shapes. This assumption allows equation ( 5 ) to be 
replaced by a set of equations of modal motion which can be solved for 
the displacement due to each mode. The bending-moment transfer function 
is then calculated by two methods: first, the mode -displacement method, 

in which the bending moment is determined directly from the curvature 
associated with the flexible mode-displacement response and, second, 
the force-summation method, in which the motion and aerodynamic loadings 
associated with the modal motions, rigid and flexible, are obtained first 
and then integrated with the applied gust loading to obtain the bending- 
moment response. Brief derivations of these modal analyses are given in 
appendix B for one- and two-dimensional turbulence. 


Matrix or Segmented-Wing Method 

Equation ( 5 ) is integrated twice to obtain an integral equation with 
z" as the unknown. In order to solve this integral equation, the wing 
is divided into n - 1 segments; thus, the equation can be written as 
an nth-order matrix equation. 

Integrating matrices of the type given in reference 6 are used to 
solve this matrix equation directly for z" and, thence, the bending- 
moment transfer function from EIz" . A brief derivation of the matrix 
method is presented in appendix C for one- and two-dimensional turbulence 


INPUT SPECTRA 


The input spectra required by equations (l) and ( 3 ) may be approxi- 
mated by several analytic expressions. (See ref. 7-) The expressions 
used herein are: 


For one -dimensional turbulence. 


- f 


1 + 3 S ^^ + a^ + 6s^a^k^ - 3 S ^^ 

1 + s^k^) (a^ + s^k^ ^ 


(3 


w 


(6) 
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For two-dimensional turbulence, 


» u ( k ,e> - ffif 


(l + 3s g k 2 )e 


(l + s 2 k 2 ) 


3/2 


K l( e 'J 1 + s 2 k 2 ) - - — 9 K 0 ^l + s 2 k 2 ^ 


2 a 2 0 2 


( a g 4- ?B gk g ) U a 2 + s 2 k 2-\ 

(.2 + s 2 k 2) 2 ° l ' 

(a^ + 6a^s^k^ -■ 3s^k^)0 - (a^ + s £ 'k^)a i+ 0^ 


(a 2 + s 2 k 2 ) 


5/2 


K ll e ^a 2 + s 2 k 2 


(7) 


where 

s dimensionless scale of turbulence, 

5/2 

L longitudinal scale of turbulence 

a arbitrary constant used to select frequency at which spectrum 

begins to deviate from an inverse second -power variation 
with frequency and begins to c.ecrease more rapidly 

On the basis of the measurements reported in reference 7, a value for 
a of 50 appears appropriate and is used in the present calculations. 
(Note that eq. ( 7 ) reduces to eq. ( 6 ) as 0 ->0.) 


DETERMINATION OF OUTPUT SPECTRA 


For one -dimensional turbulence , the oulput spectrum is obtained 
simply and directly from equation (l) for all methods considered. For 
two-dimensional turbulence , however, it should be noted from equation ( 3 ) 
that a double integration of the product of the input spectrum and the 
convolution of the transfer function is nece ssary before the output 
spectrum can be evaluated. For the exact ard. matrix methods these inte- 
grations were performed numerically because of the complexity of the 
integrands. For the modal methods, however, equation ( 3 ) is simplified 
by the separation of the frequency and space variables in the transfer 
function; that is. 
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r-1 


H(k,T)) = £ H^kbnU) 


( 9 ) 


n=0 


where r indicates the number of modes , including rigid-body modes, in 
the analysis. This separation of variables reduced equation (3) to the 
following form: 


j i 


- X I Re [M k ) H n*( k >]*w im ( k ) 


n=>0 m=0 

where for mode -displacement method 

j = r - 1 

and for force-summation method 

j = r 

and 


( 10 ) 


o w (k) = / r(e)<Jv(k,e)de 

Jq 

(11) 

r b/2 

r(e) = / 7 m (Ti)7 n (Th-6)dT) 

(12) 

J -b/2 



The force-summation method contains more terms than does the mode- 
displacement method due to the externally applied gust load. (See 
appendix B.) It should also be noted that the transfer function H n (k) 
is not identical in the two methods; for example, see equations (Bllb) 
and (B 15 b). 

The integrand of equation (l 2 ) is a relatively simple function and 
in most cases is integrated in closed form. However, the integration 
indicated in equation (ll) is performed numerically because of the com- 
plexity of the input spectrum. 


SCOPE OF CALCULATIONS 


Output spectra and mean-square values of the bending-moment response 
at the wing root and the 0. 75“ sem i s P an station for one- and two-dimensional 
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turbulence have been computed for three configurations of the simplified 
airplane (fig. l) . These configurations, defined by the variable p, 
the ratio of fuselage mass to total mass, are: 

p = 0 (all mass uniformly distributed across the wing span) 

p = 0.5 (one -half the total mass concentrated in the fuselage and 
one-half uniformly distributed across the span) 

p = 0.8 (80 percent of the mass concentrated in the fuselage and 

20 percent distributed uniformly across the span) 

In all computations, the aspect ratio was assumed to be 12.0, and 

the dimensionless scale of chordwise turbulence s = -t- was assumed to 

c/2 

be 100. In order to facilitate the computations, the reduced frequency 
associated with the first free-free symmetrical mode of the dynamic 
system involving wing bending and the wing aerodynamic mass parameter 

were assumed to be invariant with n . These assumptions reduced the 
computational labor appreciably but resulted in each of the three air- 
plane configurations having a different total airplane mass. Since, 
however, the purpose of this study was to evaluate methods rather than 
to show trends, this difference was not considered a serious drawback. 

The values chosen for the constants, k^ = 0.2 and ^ = 100, cor- 
respond to those of an aircraft which has a high degree of wing flexi- 
bility and low aerodynamic damping. Such a system, it was felt, would 
provide a severe test of the adequacy of the approximate methods to be 
evaluated. Since the wing mass parameter was held constant, the airplane 
aerodynamic mass parameter varied with p. - that is, for \i = 0, k = 100; 
for [i = 0.5; k = 200; and for = 0.8, < = 500 . 

Aerodynamic strip theory was used in akL of the studies since a more 
refined theory would have further complicated the already complex solu- 
tions and probably would not have Influence! appreciably the conclusions 
regarding the relative accuracy of the various methods. Also, the 
unsteady lifts associated with the gusts anl the resulting wing motions 
were included in the analyses. 

All the results presented herein are ii nondimens ional form v The 
output bending -moment spectra and the mean- square results have been 
divided by the square of the bending moment due to a sharp -edge gust of 
unit vertical velocity which is constant across the span. 

The modal methods require a finite numoer of modes to approximate 
the airplane motions. For the mode -displacement method, rigid-body 
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translation and the first natural wing bending mode were used. For the 
force -'Summation method, the following three approximations were 
investigated: 

(1) Rigid-body translation and the first natural wing bending mode 
(same as for the mode -displacement method) 

(2) Rigid-body translation and a parabolic approximation to the 
first free -free wing bending mode 

(3) Same as approximation (2) except that the frequency associated 
with the natural wing bending mode was used instead of the approximate 
Rayleigh frequency 

The matrix method requires that the wing be divided into a number 
of equally spaced stations across the span in order to perform the 
matrix integrations. For the present study, five such stations across 
the wing semispan were used. 

For purposes of comparison, bending -moment results were also calcu- 
lated for the simplified airplane considered to be rigid. For these 
calculations the force-summation method was used with the airplane motions 
limited to rigid-body translation only. A brief derivation of this pro- 
cedure is given in appendix D. 


ACCURACY OF NUMERICAL INTEGRATIONS 


Inasmuch as the integrations involved in calculating the two- 
dimensional output spectra were performed numerically (see the section 
entitled "Determination of Output Spectra"), a check on the numerical 
accuracy was deemed necessary. 


Exact Method 

For the exact method a suitable check on the numerical accuracy of 
the two-dimensional expression was obtained by calculating the one- 
dimensional output spectra from the two-dimensional expression (see 
eq. (3)) and comparing the results with those from the one -dimensional 
expression (eq. (l)). This one -dimensional check calculation is obtained 
by substituting $ w (k) for $v(k,9) in equation (3). Thus, any dis- 
crepancies between the results of the two procedures can be traced 
directly to the numerical integrations of the two-dimensional expression. 
A comparison of these results is given in figure 2. 
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From part (a) of figure 2 it is evident that both methods give 
nearly the same results for the spectra of root bending moment. It can, 
therefore, be expected that the accuracy of the root -bending -moment 
results for the two-dimensional gust structure are satisfactory. 

From part (b) of figure 2 it is evident that both methods do not 
give the same results for the spectra of bending moment at the 0.75- 
semispan station. It can, therefore, be expected that the results at 
the 0.75-semispan station for the two-dimens Lonal gust structure will 
not be very accurate. The inaccuracy of these results can, in large 
part, be attributed to the fact that the 0. 7 :> -semispan station is located 
between two of the equally spaced stations used in the numerical integra- 
tion procedures. This location, unfortunate Ly, places the numerous dis- 
continuities in the slope of the transfer functions (see fig, 5) to be 
integrated midway between the stations. At the time of this discovery 
it was considered impractical to change to another spanwise station; 
thus, the results at the 0.75“ se rcispan station calculated by the exact 
method will not be presented for the two-dimensional turbulence condition. 


Modal Methods 

For the modal methods the integrations to be performed in order to 
evaluate the output spectra are defined by equations ( ll) and ( 12) . The 
first integration, equation (12), is of a form which can be integrated 
analytically or accurately by numerical procedures. Thus, the subsequent 
integration, equation (ll), which is performed numerically, should also 
be accurate. 


Matrix Method 

The numerical integrations used to calcilate output spectra for the 
matrix method are the same as those used in the exact method. Figure 5 
shows that excellent agreement exists between the exact and matrix trans- 
fer functions; thus, the numerical inaccuracies in this method are pre- 
dominantly the same as in the exact approach. For this reason, the 
matrix results at the 0.75- Ber ni s P an station Tor two-dimensional turbulence 
are not presented. 


EVALUATION OF APPROXIMATE PROCEDURES 


In this section the various approximate procedures are evaluated 
from the output-spectra and mean-square bending -moment results for one- 
and two-dimensional turbulence. In order to facilitate the discussion, 
the approximate procedures are divided into vwo groups: first, the more 
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accurate methods consisting of the force -summation method based on a 
natural mode and the matrix method; second, several simplified modal 
methods consisting of the mode -displacement method based on a natural 
mode and the force -summation method based on a parabolic mode with the 
corresponding Rayleigh-Ritz frequency and with the natural frequency. 

For all solutions the results for the three model configurations 
studied show the same type of comparison; thus, the output spectra are 
presented for only one representative configuration, namely, p = 0.5* 


Force-Summation and Matrix Methods 

Statistical bending -moment results obtained by the force-summation 
method using one natural mode and the five-station matrix method are 
evaluated by comparison with the exact-method results. In order to per- 
mit evaluation of the influence of the second symmetrical bending -moment 
resonance peak on the results, the output spectra have been calculated 
by the exact and the matrix methods over a larger range of reduced fre- 
quencies than is usually considered necessary. The output spectrum for 
the airplane assumed to be rigid is also presented for comparison. 

Bending moment at the wing root . - Spectra of the bending moment at 
the wing root for a one -dimensional gust structure are presented in 
figure 4(a) . It is readily apparent that both of the approximate methods 
yield spectra which are almost the same as those for the exact-method 
results. Also, the results for the rigid airplane agree very well with 
the results for the flexible airplane calculated by the exact method for 
small values of the reduced frequency where dynamic flexibility effects 
are unimportant. 

Spectra of the bending moment at the wing root are compared in 
figure 4(b) for the case of a two-dimensional gust structure. The agree- 
ment is again seen to be very good throughout the frequency range. 

For both the one -dimensional and the two-dimensional gust structures, 
the peak of second symmetrical bending mode, located at k = 1.12, is so 
small that it is barely evident on the plots and, thus, has little effect 
on mean-square results. 

Bending moment at 0.75-semispan station .- Spectra of the bending 
moment at the 0.75-semispan station for a one -dimensional gust structure 
sire compared in figure 4(c). The two approximate solutions agree very 
well with the exact -method results. The force -summation solution, of 
course, does not give the second-mode resonant peak, and the matrix 
method gives it displaced to a somewhat higher frequency. In this high 
frequency range there is no experimental confirmation of the shape of 
the input gust spectrum; thus, the significance of the second-mode peak 
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is questionable. For mean-square results th: s peak can be neglected, 
but for higher statistical moments, such as are needed to predict the 
time expected to exceed a given b ending -momer.t level, these secondary 
peaks may be important. 

For the two-dimensional case it was found that the exact- and 
matrix-method results for the 0.75-semispan station were inaccurate due 
to numerical integrations. (See section entjtled "Accuracy of Numerical 
Integrations.") Hence, figure 4(d) presents only the force-summation 
and rigid-airplane results. Since the numerical integration errors pres- 
ent in the exact and matrix methods are not present in the force -summation 
method, there is no reason to believe that this method is any less accu- 
rate for this case than for any other case considered. Figure 4(d) lends 
support to this statement, by showing good agreement between the force- 
summation method and the rigid analysis at lew frequencies. It should 
be remembered that the numerical inaccuracies in the exact and matrix 
methods occurred in the determination of the output spectra and not in 
the transfer -function calculations. Thus, as was shown in figure 5, the 
matrix method yields excellent two-dimensional transfer-function results. 


Simplified Modal Methcds 

Statistical bending -moment results obtained by simplified modal 
approaches are evaluated by comparison with the force -summation results 
based on a natural mode. (This method has been shown to be accurate in 
the preceding section.) The simplified modal methods consist of the 
mode -displacement method based on a natural node and the force -summation 
method based on an approximate parabolic mode . Two different procedures 
are used in the latter method: the Rayleigh-Ritz procedure in which the 

Rayleigh frequency associated with the appro* imate mode is used, and the 
modified Rayleigh-Ritz procedure in which the frequency associated with 
the natural mode is used. 

Bending moment at the wing root .- Spectra of the bending moment at 
the wing root calculated by the modal methods for a one -dimensional gust 
structure are shown in figure 5(a) . Compared with the force -summation 
method based on a natural mode, the mode -displacement method yields inac- 
curate results in the low frequency range. 1 he force -summation Rayleigh- 
Ritz approach based on the parabolic mode results in a shift of the res- 
onance peak, as would be expected, because of the Rayleigh frequency. 
Although, this approach yields a fairly good result for the mean-square 
bending moment, it would be a very poor choice if higher statistical 
moments were desired. The force -summation modified Rayleigh-Ritz approach 
based on the parabolic mode is almost as accurate as the results when 
the natural mode is used. 

The results for the two-dimensional gust structure are presented in 
figure 5(b). The same trends can be seen as for the one-dimensional case. 
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Bending moment at the 0 . 75 -semispan station ,- Spectra of the bending 
moment at the 0 . 75 -semispan station obtained by the modal methods are 
shown in figures 5 (c) and 5 (d) for the one- and two-dimensional gust 
structure , respectively. It is again apparent that the mode -displacement 
method is very inaccurate in the low frequency range, whereas the force- 
summation modified Rayleigh-Ritz approach is very accurate throughout 
the frequency range shown. Results were not calculated by the force- 
summation Rayleigh-Ritz approach based on the parabolic mode for this 
station. 


EFFECTS OF SPANWISE VARIATIONS OF TURBULENCE 


In this section the effect of the spanwise variation of vertical 
turbulence on the bending -moment response is assessed. First, the effects 
on the output spectra are shown in figure 6. Part (a) of this figure 
shows the spectra at the wing root as determined by the exact method. 

For the configurations with p = 0-5 and 0.8, it is apparent that the 
main effect of the spanwise variation of turbulence is an appreciable 
reduction in the bending-moment response as compared with the one- 
dimensional case, particularly in the vicinity of the fundamental mode 
resonant peak. For the configuration with p = 0, however, the two- 
dimensional turbulence results in an increase in root -bending -moment 
response because the response produced by the one -dimensional turbulence 
is zero. Part (b) of figure 6 shows the spectra at the 0. 75- sem l s P an 
station as determined by the force-summation method. (The exact-method 
results have been shown to be inaccurate for this station because of 
numerical integrations.) The same results are apparent at this station 
as at the root station. 

Second, the effects of spanwise variations of turbulence on the 
mean-square bending -moment results are shown in figure 7- As would be 
expected, the same effects shown on the spectra are further amplified 
in the mean-square results; namely, a decreased response for p = 0.5 
and 0.8 and an increased response for p = 0. These results are in 
agreement with those given in reference k for root -bending -moment 
response. 

Finally, the effects of spanwise variations of turbulence on an 
amplification factor, defined as the ratio of the mean-square bending 
moment of the flexible airplane to the mean-square bending moment of the 
rigid airplane, are shown in figure 8* Again the same general conclu- 
sions that were reached for the results shown in figures 6 and 7 are 
evident. In addition, it should be noted in the low p range (p -» 0) 
that the trend of the amplification factors for the one-dimensional gust 
is not the same as that for the more realistic two-dimensional gust. 
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CONCLUSIONS 


A study has been made of the relative accuracy of several approximate 
procedures for calculating the symmetrical bending -moment response of 
flexible airplanes to continuous isotropic turbulence. These approximate 
procedures have been applied to a simplified airplane which consists of 
a uniform beam with a concentrated fuselage mass at the center. The 
conclusions drawn from this study are: 

1. The force-summation method based on one natural bending mode gives 
very good results compared with the exact soLution. 

2. The matrix method based on five stations across the semispan gives 
very good results compared with the exact soLution. 

3. The mode -displacement method based on a natural mode yields inac- 
curate results; however, this relatively simple method can be useful in 
trend studies involving variations in wing flexibility. 

4. The force -summation method, based on an approximate parabolic 
mode and the Rayleigh frequency, loses little accuracy in mean-square 
results. However, the approximated natural frequency causes a shift in 
the fundamental mode resonant peak; conseque itly, higher statistical 
moments will be in appreciable error. 

5. The force -summation method based on in approximate parabolic mode 
loses little accuracy if the natural frequency is known. 

In addition to the evaluation of approximate methods, the effects 
of spanwise variations of vertical gust velocity were also studied. From 
this study it was found that the inclusion of spanwise variations of tur- 
bulence results in a decreased response if most of the airplane mass is 
in the fuselage and in an increased response if most of the mass is in 
the wing. Also, the exclusion of spanwise variations of turbulence in 
trend studies, such as the trend of amplification factor with mass ratio, 
may lead to erroneous conclusions. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., October 2k, 19o8. 
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APPENDIX A 

EXACT ANALYSIS FOR SIMPLIFIED MODEL SUBJECTED TO 
SYMMETRICAL PART OF RANDOM-GUST DISTURBANCES 


The differential equation defining the symmetrical motion of one- 
half the simplified elastic airplane consisting of a uniform beam with 
a concentrated mass at its center and subjected to a concentrated sinus- 
oidally varying gust at station t) (see fig. l) is 


EIz”" = -m(y)z - 1 mp,z6(y-0) - C(k)c 7 qc + 0(k)c 7 qc e^&Cy-q) 

2 r ‘•a U ‘•a, U 

(Al) 

Equation (Al) is a specific form of equation ( 5 ) . Since the applied 
force is varying sinusoidally, the steady-state response also varies 
sinusoidally with the same frequency and thus a solution of the form 

z(y,t) = z(y,a))e in)t can be postulated. Nondimensionalizing equation (Al) 
according to the definitions given in the symbol list yields 


= j£_|i + R 8 (y-0)] z - 2i z + £ 6(y-rj) (A2) 




k a 2|t w U 


By using the root boundary conditions z'(o) = z m (o) = 0, the solu- 
tion of equation (A2) is obtained by the Laplace transform as 


Z = 


_ z(0) 


cosh ay + cos ay (sinh ay - sin ay) + 


) 


Z n (0) 

2a 2 


(cosh ay - cos ay) + — — [sinh a(y - tj) 


(A3) 


where 


sin a (y - Ti)]7(y-^i) 
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when the second and third derivatives of equation (A3) are taken and the 
tip boun dar y conditions, z"(l) = z'" (l) =0 are used, the unknowns 
Z(0) and z"(0) are determined. Substituting these root conditions 
into equation (A3) and taking the second derivative yield 



where 

= cosh a + cos a 
C 2 = cosh a - cos a 
= sinh a + sin a 
= sinh a - sin a 
F-^ = cosh ay - cos ay 
F 2 = sinh ay + sin ay 
Fj = cosh ay + cos ay 

D = 2 [A(l + cosh a cos a) - a^(sinh a cos a + sin a cosh a)J 

Since a = g - ih is a complex number, further manipulations are required 
to put the expression for z" into a form suitable for computation. The 
final form of z" is 
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z"(M) = H z «^ c [5 R (k,Ti) + (A5) 

where 

g 000 

2(Qi 2 + Q 2 2 ) k a 2ft w 

%( k > 4 ) = Q9P57 + ^10^68 “ ^11^15 ” Qi2%£ + a 3^57 + hP 68^°0 
HjCk^Tj) = QxqP^y “ ^68 ” ^12^13 + ^11^2^ + A 3 ( gP 57 “ hP 68)-/( c 0 

and 

Q9 = + QgQg 

^10 = ^ 2^5 _ q i q 6 

Q ll = ^1^7 + Q 2^8 
Q]_2 ~ Q2Q7 “ Q1Q8 

Q 5 = d T P i3 “ d 8? 24 + d 9 P 57 " d lO P 68 + d ll p 13 “ d l2 P ^4 

06 = dyPgij. + + d 9P68 + d 10^57 + d ll%q: + d i^i3 

Q7 = a^— - d 2 P 24 + d 3 P 57 “ d 4 p 68 + d 5 p l3 “ d 6 P 2£ 

% ~ d i p 24 + d 2 p ^3 + d 3 p 68 + d 4 p 57 + d 5 P 24 + d 6 p i3 

d i " - p 1 3 (s 2 - h2 ) + 2 e hP 2it 

->2 - -F^g 2 - h 2 ) - 2ghP 15 
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d 5 - PjtU 2 - h2 ) - 2ghP 2 4 - / — j (e%f h hp 68) 

■V) - p 24(® 2 - h2 ) + 2ghF y - ■ (" g 2 2 h 2 ^ (^68 • hp 57 ) 
d 7 = -^(g? - h 2 ) + 2 ghP 5g 

d 8 - -^fe 2 - h2 ) - 2ghP 57 

<9 -(yh?) (**♦“*) 

dl ° “ (?T1?) ' “” 57 ) 

d ll = P 57 (s 2 - h 2 ) - 2ghP 68 - d 5 

di2 = p 68( g2 " h2 ) + 2ghP 57 " dl + 

2 (Ql 2 + Q 2 2 ) 

A 3 2 2 — 

g 2 + h 2 

Ql = Mi - g(g 2 - 3 h 2 )q 2 + h(h 2 - 3 g 2 )q 3 
02 = Aq 4 + g(g 2 - 3h 2 )q 5 + h(h 2 - 3 g 2 )q 2 

qi - 1 + P 1 P 3 + P 2 P 4 
q 2 = p 3 p 5 + P 4 P 6 + P 1 P 7 ■ P 2 P 8 
■J 3 - P 4 P 5 - P 3 P 6 - P 2 P 7 - P 1 P 8 
84 - p 2 p 3 - P 1 P 4 



23 


= cosh g cos h 
P 2 = sinh g sin h 
= cos g cosh h 
P4 = sin g sinh h 
P^ = sinh g cos h 
P^ = cosh g sin h 
P^ = sin g cosh h 
Pg = cos g sinh h 


P^ = cosh gy cos hy 
P 2 - sinh gy sin hy 
P^ = cos gy cosh hy 
P4 = sin gy sinh hy 
P^ = sinh gy cos hy 
P^ = cosh gy sin hy 
Py = sin gy cosh hy 
Pq = cos gy sinh hy 


= cosh g7 cos h7 
P 2 = sinh g7 sin h7 
P^ = cos g7 cosh h7 
P4 = sin g7 sinh h7 
P^ - sinh g7 cos h7 
P^ = cosh g7 sin h7 
P rj = sin g7 cosh h7 

A 

Pg = cos g7 sinh h7 


P-j_ = cosh ga cos ha 

\ 

P 2 = sinh ga sin hx 
P^ = cos ga cosh ha 
P4 = sin ga sinh ha 
P^ = sinh ga cos ha 
P^ = cosh ga sin ha 
Py = sin ga cosh ha 
Pg = cos ga sinh ha 


^mn 


p m + 

p n 

^rnh 

= 

P m - 

P n 

7 = 

1 

- n 


a = 

y 

- tJ 



The bending -moment response is obtained from 



Nondimensionalizing equation (a 6) by the bending -moment response due to 
a sharp-edge gust yields 


M 


B 


m B,s 


0(k) 

^(Q] 2 + Q2 2 )^cp 


H R (k,fi) + iH- 



= H R (k,fj) + iH I (k,i]) 


(A7) 


The bending -moment transfer -function response to a one -dimensional 
gust field may be obtained from the two-dimensional transfer function as 
fo Llows: 


H(k) 



H(k,Tj)drj 


(A8) 
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Integrating equation (a 8) yields 

H(k) = Hpj(k) + !%(}:) 


(A9) 


where 


H R (k) 


v g 2 + h 2 / 


Q 9 P 5T + ^lO 1 ^ ~ °'11 P 13 " Q 12 P ^4 + 


A j( gP 57 + **68) 


0(k) 

/ 2 p\ 

H Q i + Qg y y cp 


and 


H T (k) = 


lg 2 + h 2 / L 


q 10 P 57 ” q 9 P 68 " °"12 P 13 + Q ll P ?r + 



P 13 gP 57 + **68 


P — 

24 


_ gP 68 " ^57 



gP l3 + 


hP 


24 


p 68 = gP 24 - ^ 
p 57 = gp i3 + ^24 
P 68 = gP 24 ' **13 


All other terms are as defined after equation (A5) for the two-dimensional 
case . 
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APPENDIX B 


MODAL ANALYSES FOR SIMPLIFIED MODEL SUBJECTED TO 
SYMMETRICAL PART OF RANDOM-GUST DISTURBANCES 


In this appendix, solutions of the differential equations defining 
the symmetrical motion of the simplified model illustrated in figure 1 
are obtained by two modal techniques, namely, the mode-displacement 
method and the force-summation method (ref. 5)* Within the framework 
of these two methods, solutions are obtained by the Rayleigh-Ritz pro- 
cedure based on an exact free-free bending mode, and, in addition, for 
the force-summation method, on an approximate parabolic mode. 


Displacement Analysis 

The determination of modal deflections is a common initial step in 
both the mode -displacement and the force -summation methods. The deflec- 
tion expressions are a function of mode shape and thus give different 
results depending on whether exact- or approximate -mode shapes are 
employed . 

The Rayleigh-Ritz equations defining the displacement in the assumed 
modes of an unrestrained elastic airplane subjected to an arbitrary 
external loading are given in reference 5 . For the simplified airplane 
of the present study restricted to vertical rigid-body displacements and 
deflection in an assumed fundamental free-free symmetrical bending mode, 
these equations can be written in the terminology of the present report 
as 



(B2) 
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where 

z 0 

Z 1 

W 


generalized 

generalized 

symmetrical 


coordinate for rigid-body translation 


coordinate for first free-free wing bending mode 


part of complex amplitude of gust velocity 


In writing these equations, the vertical displacement at any point on 
the wing was expanded in terms of the rigid-body vertical-motion mode 
and some mode ^(y), representative of the fundamental free-free mode 

of the simplified airplane, as 


z = z 0 + z 1 | 1 (y) 


(B3) 


It was further assumed that the mode |]_(y), whether exact or approximate, 
satisfies the orthogonality relation with the rigid-body mode, or that 

j~,b/2 

/ n4,(y)li(y)<iy = o 

j -b/2 

If the second term of equation (B2) is multiplied and divided by the 

P 13 / 2 r g 

expression J ni^Cy) [ |^(y)j dy, it can be written as 


Zn < 


-b/2 

n b^2 q 

/_ b/g EI [ Sl " ( y)l dy 
/-b/2 "t (y) [*l (y) ] ^ 


J m t (y)[li(y)] 2 dy 


If the square of the approximate Rayleigh frequency o>^ a is substituted 
for the term in braces, this expression becomes 

.b /2 


u/ ^ 2 

m t (y) [li(y)J dy 


Using this expression and nondimen sionaliz in 5 equations (Bl) and (B2) 
give 
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2iCOO/- 
*w k V 


K ^l) 


«„k 2 U 


(B4) 



m(y) 


^t(y) 

m t (0) 


Solving equations (B4) and (B5) simultaneously for zq and z^ gives 
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If an exact free-free bending mode for the system is used rather 
than an approximate shape, the Rayleigh frequency o>^ a is identical to 

the exact fundamental -bending-mode frequency for the system o>p and, 
thus, kp = k^ a . Where the exact frequency is known but the exact mode 

is unknown, an increase in accuracy can normally be obtained by using 
the exact frequency in place of u>^ a . This procedure is referred to as 

a modified Rayleigh-Ritz approach in the present report. 


Bending Moments Due to a Two-Dimensional Gust Distribution 

Mode-displacement method .- In the mode --displacement method, the 
bending moment is defined by simple beam theory as 


r-1 

M b = El(y x ) £ z n ! n "(yi) 0*9) 

n=0 

In the present study, El is a constant anc. So"(^l) '*' s e< ^ ua -*- to zero > 

thus, 

M b = EEz 1 | 1 n ^y 1 ) (BlOa) 

or 



(BlOb) 
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Dividing equation (BlOb) by the bending moment due to a sharp -edge 
gust and substituting equation (B 7 ) for zp give 


= Ho L + H X L -L 
M B, s S g ’ 


(Blla) 


where 


H, 


Hi 


2iK 4 


kC(k) 


2 

K w k Ra 

k Ha 2 y C p A 


k 2 0 . kC(k) 

k 2 f "i 
K 1 *1 1 

(yi)0OO 

^1^-Ra ^v^Ra 

-^Ra ^cp ^ J 


( Bllb ) 


If the procedure described in reference 4 is used, the spectrum of the 
output bending-moment ratio is written as 


Ks) 2 

where 


and 


H o| 2<t w 00 (k) + | H l| 2<3) w 11 ( k ) + 2Re ( H o H l*) 5> w 01 ( k ) 




7 m (n)7 n (ti+0)dTi 


7 0 ( H) 


= 1 

= it ( tj) 


(B12) 


(B13) 


Equation (B13) is a slightly modified version of equation (ll) . 

Equation (B13) can be used with either an exact bending mode or an 
approximate mode which satisfies the orthogonality condition with the 
rigid-body mode. Unless the second derivative of the approximate mode 
is a very good approximation of the second derivative of the exact mode, 
however, the bending -moment result cannot be expected to be very accurate. 


Force-summation method .- In the force- summation method, the inertia 
and the aerodynamic loading due to wing motions in the various modes are 
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integrated together with the external gust loading to obtain the desired 
bending moment as follows: 


M 


B 



l m + l g)(y - yi) d y 


(B14) 


where l ^ and Z m are the inertia and aerodynamic loadings associated 

with the deflections of the rigid body and fundamental modes and l „ is 

& 

the symmetrical part of the loading due directly to the gusts and they 
are defined by the following equations; 

l± = m t (y)cn 2 (z 0 + z^Cy)) 

c ( k )( z o + zi6i(y))i^ 


c cq 



U 


0(k)w 


Substituting the previously derived expressions for Zq and 

(eqs. (B6) and (B7)) into equation (Bl4) and manipulating the result 
gives y after dividing by the bending moment due to the sharp-edge gust, 


where 


M- 


B 




- H 0 L g + H l L g,l + H 2 L g,2 


(B15a) 


H 0 = (Tjl + iT 2 )(T, + 1T 4 ) 

% ■ (T l + iT 2 )(T 5 + ITg) 
= A = Ty + iTg 

L g,2 = \_ (y - h) w d y 

v v-. 



(B15b) 


yi 
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and 



G = 2k[G(k) + & 

*wL 2 . 


F = 2k F(k) 
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pi 

I 2 = / y6 x (y)ay 

J y ± 

As = k 5 + 1^(0)] 2 

When the procedure given in reference 4 is used, the spectrum of output 
bending -moment ratio is written as, 

= |H n | 2 ® (k) + |H,| 2 4 (k) + |Hp| (k) + 2Re(H n H *U (k) + 

( M B s ) 2 ° °° 1 11 2 22 V 0 1 / W 01 

2Re(HoH 2 *)<tv 02 (k) + 2Re (HiH 2 *) ^(k) (Bl6) 

where the functions ^^(k) are as defined in the preceding section 

and 7 2 (n) = | (j\ ~ 9 1 + ( “ 1 ~ 9jJJ ('^“^ 1 )] • E( l uation (Bl6) ls 

an expanded version of equation (ll) and can be used with an exact bending 
mode or an approximate mode which satisfies :he orthogonality relation 
with the rigid-body mode. 


Bending Moments Due to a One-Dimensional Gust Distribution 


Mode -di splacement method . - In the mode- displacement method for the 
one-dimensional gust distribution, the integrals defining Lg and Lg^p 
can be evaluated exactly. These integrations yield 


L = 2w 

g 

L g, 1 = Zvfy > 


(B17) 


Substituting equations (B17) into equation (dll) gives 

Mb 


m b ,; 


= (2H 0 + 2K 1+ H 1 w 


(B18) 


Then the output bending -moment spectrum becomes, 


%B^ k ) 

( M B,s ) 2 


2Hr 


2K^H 1 | 2 + (B19) 
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Equations (Bl 8 ) and (B19) can be simplified somewhat by combining terms, 
but this was not done since Hq and Hq for the two-dimensional case 
had to be determined in any event. 

Force-summation method .- In the force -summation method, for the one- 
dimensional gust distribution, the integral expression for Lg^ eva l u_ 

ated exactly yields 

L g,2 = y cp w (B20) 

Substituting equations (B20) and (BI 7 ) into equation (BI 5 ) gives 

i = w( 2 H 0 + 2Kj^Hp + y cp H 2 ) 

Mb,s v ' 

The output bending -moment spectrum for the one-dimensional gust 
distribution can then be written as, 

= [4|H q |2 + 4^1^12 + y^l^l 2 + 8^(1^*) + 

4y cp Re (Hq^*) + <D w (k) (B22) 
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APPENDIX C 

MATRIX METHOD FOR SIMPLIFIED MODELS SUBJECTED TO THE 
SYMMETRICAL PART OF A RANDOM-GUST DISTURBANCE 


The differential equation defining the symmetrical motions of one- 
half the simplified wing-fuselage combinaticn treated in this paper 
(fig. l) and subjected to a sinusoidal gust velocity at station T| on 
the span is given by equation (A2) as 

z"" = *L[l + R6(y-0)] z - 2i — ^ z + £ 8(y-tj) (Cl) 

a w K a *w 

Integrating equation (Cl) twice by using twc of the four appropriate 
boundary conditions; namely z" = z MT =0 £t y = 1 gives 

z" f_ J_ |k 2 [l + BB(y-O)] z - 2ik z + 3^1 § 6(y-rj)|dy dy 

a y y i J 

(C2) 

By definition 


■jr* 


dy dy 


Using the boundary condition that the slope of the wing bending deflec- 
tion at the fuselage center line (y = 0) must be zero for symmetrical 
motions in this definition yields 


z = z 0 + f f 2" dy dy 

J oJ o 


(C5) 


Substituting equation (C3) into equation (C2) gives 
z” = U- [ f (L?[l + RB(y-O)] - 2ik ^l(z 0 + ^ T* 2" dy dy^ 

k a J y J y V *w \ jqJ o y 


£ 6(y-fj)]dy dy 


*w U 


(C4) 
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Dividing the beam into n - 1 segments, equation (C4) can be written as 

1 nlrl 


an 


nth-order matrix equation. Designating \J-l\ = / > = 

J v 



y^ y 


[I 2 ] = J' > an d [II 2 ] = J J as integrating matrices for the contin- 
uous functions in equation (CP) and integrating the delta functions 
exactly allows equation (Cl+) to be written as 


{*"} = %“ l ] 


"a 


1 - 2i 


. C(k) 


kK w 


‘o| + p 1 ^] f + ^r~ n{ (f i - y)/^-y)} 

j i k a K w 


(C5) 


Under the assumption of symmetrical bending deflections, the shear at 
the fuselage center line must be zero; thus, 


, r h / 2 

(EIz 'Vo = °=/ 0 ( 


H + l m + l g) d y 


(c6) 


Equation (c6) can be written in matrix notation as follows: 


0 = 


Li] 


i - 21 SiH 

k *w 


+ R 


L x sl 


|{z 0 } + [lI 2 ]{f"}} + ^> i (C7) 


where |_lj is a row matrix of integrating factors for the range 0 to 1 
and |IjJ = [l 0 0 0 Oj for n = 5* 


Solving equation (C 7 ) for |^o| gives 


N - i 


[I] 


1-21 


1 C(k) 
k *V 


+ R 


P 5 ] 


[II 2 ] {2") + g{l} (c8) 


where 


Ki = - 


LiJ 


1-21 


. C(k) 


kt V 


R 


M 




and [i] and are composed of identical rows of [ij and of 

|JEjJ , respectively. 
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Substituting |zq j from equation (c8) 
after rearranging terms , 


:.nto equation (C5) yields. 



From equation (C9), z M is obtained by inverting the matrix premul- 
tiplying {£"} at any or all of the n stations due to a sinusoidal 
gust of intensity w applied at any station tj on the wing. Note that 
T] need not be at the stations into which the wing is divided; it may be 
anywhere on the wing. In the computations performed in the present 
report, only 5 stations were used but 10 equally spaced values of r\ 
were required to obtain reasonable accuracy :n subsequent numerical 
integrations . 


In the computations, the following integrating matrices from refer- 
ence 6 were used: 


LiJ = 

L 1 

4 

2 

4 

8 


2 

2 k 

44 

88 

54 


0 

2 

21 

60 

25 

m - (0 2 f 

0 

0 

0 

52 

16 


0 

0 

-1 

6 

7 


0 

0 

0 

0 

0 


~ 0 

0 

0 

0 

o" 


7 

6 

-1 

0 

0 

m - (0 2 f 2 

l6 

52 

0 

0 

0 


25 

60 

21 

2 

0 


54 

88 

44 

24 

2 
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The bending moment is obtained from 

El 


M b = EIz" = 


/ b\ 

\ 2 / 


d) z - 

2 


(CIO) 


Nondimensionalizing equation (CIO) by the bending -moment response due to 
a sharp -edge gust yields 


= a — z" = H B (k, rj) + iHjUjij) 

Mb.s 2y„ — 


(Cll) 


The bending-moment transfer-function response for a one -dimensional 
gust field may be obtained by simple numerical integration over t) of 
the two-dimensional transfer function; that is, 

H(k) = f H(k, f})drj (C12) 

J 0 
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APPENDIX D 

ANALYSIS FOR SIMPLIFIED RIGID MODEL SUBJECTED TO THE 
SYMMETRICAL PART OF A RANDOM-GUST DISTURBANCE 


In this appendix, an exact solution is made for the symmetrical 
response of a simplified rigid airplane to random vertical gust excita- 
tion. (See fig. 1 .) The basic differential equation defining vertical 
motions of this system is 


0 = l i + l m + l g ( D1 ) 

(which is a specific form of eq. (5)) where l ± is the inertia loading 
per unit span and Zg and Z m are, respectively, the aerodynamic 

loadings per unit span due directly to the symmetrical part of the gust 
and those due to wing vertical motions. Since the flexibility term is 
missing in equation (Dl) , it can be solved exactly for the system motions. 


Displacement Analysis 


Under the assumptions of the present study (see appendix B) , equa- 
tion (Dl) can be rewritten as 


/ \ n by/2 

mt(y)z 0 + C L J3q i Q _J c^cq g 0(k)dy = 0 


(D2) 


If sinusoidal excitation is assumed, equation (D2) is solved directly 
for the rigid-body displacement. This has been done in reference 3, and 
in the terminology of the present analysis the equation is written in 
nondimensional form as 


2 0 


1 

0(k) 

k 2 

-« + 2i C(k > 


k 



(D3) 


Bending -Moment Analysis 

The bending moment is obtained from the following simple integration 
of the loadings (force -summation method): 
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where, for the rigid airplane, 

l i = § mJy^ZQ 
2 

cj qin> C(k) 

~t a 2 — 

m U z 0 

! g = o Ia qcf*(k) S 


Substituting equation (D3) into equation (D4) and dividing by the bending 
moment due to a sharp-edge gust results in the following expression: 


Mr 

— s- = HfX + HoL ct 0 
M B,s 2 S>2 


(D5) 


where 


H 0 = 



1 


000 

yep 


kKn[k(K + 1) + 2G(k) - 2iF(k)] 
[k( k + 1) + 2G(k)] 2 + l*[F(k)] 2 


If the same procedure employed in the analysis for the flexible airplane 
is followed in the analysis for the symmetrical response of the rigid 
airplane, the spectrum of the output bending -moment ratio for the two- 
dimensional gust distribution is 


hm(fc>] r 

(>M 2 


l H °| 2 *“00 (k> + l H 2 | 2 % 2 <k) + 2Re ( H 0 H 2*)«v, 20 (K) 


(D6) 


For the one -dimensional gust distribution - that is, 
the span - equation (D6) reduces to 


r 

("»,.) 2 


1 + 2H 0 | 2 $ w (k) 


w constant along 


(D7) 



40 


An exact solution can be obtained for the output bending-moment 
spectrum of a rigid airplane traversing a random two-dimensional gust 
field with the distribution characteristics defined by equation (7) . 
solution, however, is highly subject to small difference difficulties 
and since the functions involved have not been tabulated with a suffi 
cient degree of accuracy, no results have been computed which may be 
used as a check comparison on the numerical integrations performed • 


The 
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H- 

(a) y, = 0; exact method. (*) 



H- 

(*) y, = 0-75; force -summation method. 

Figure 8.- Comparison of mean-square bending-moment ratio for one- and 

two-dimensional turbulence. 
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